Thermalization of an impurity cloud in a Bose-Einstein condensate. 
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We study the thermalization dynamics of an impurity cloud inside a Bose-Einstein condensate at 
finite temperature, introducing a suitable Boltzmann equation. Some values of the temperature and 
of the initial impurity energy are considered. We find that, below the handau critical velocity, the 
macroscopic population of the initial impurity state reduces its depletion rate. For sufficiently high 
velocities the opposite effect occurs. For appropriate parameters the collisions cool the condensate. 
The maximum cooling per impurity atom is obtained with multiple collisions. 



In a Bose-Einstein condensate (BEC) many atoms pop- 
ulate the same quantum state. This system is, therefore, 
suitable to realize an atom laser That is obtained 
creating atoms in an untrapped or antitrapped level by 
means of a rf or Raman outcoupling. Elastic collisions 
between impurity atoms and the parent BEC have been 
observed at MIT 0. These collisions can deplete || 
and heat the outcoupled atoms. In this work we study 
the thermalization dynamics of the impurity atoms. We 
find the following main results: i) we evaluate the impu- 
rity thermalization and depletion time for some values of 
the temperature and of the initial impurity energy. The 
thermalization is not suppressed below the Landau crit- 
ical velocity g for a finite condensate temperature T, 
according to a previous work ii) below the Landau 
critical velocity the macroscopic population of the atom 
laser state slows down the impurity depletion, whereas, 
for sufficiently high velocities the macroscopic popula- 
tion enhances the scattering rate, as observed in Ref. |^; 
iii) the collisions cool the condensate if the impurity ini- 
tial energy ej„ is below its thermal equilibrium energy, or 
even if this condition is not fulfilled, provided the thermal 
energy is sufficiently smaller than the chemical potential. 
These effects can be observed creating a succession of 
atom laser packets and observing a reduction of the im- 
purity depletion rate in the last packets. The succession 
of packets can be used to reduce the condensate temper- 
ature, as proposed in Ref. [|. In Ref. || we evaluated the 
scattering rate of an impurity inside a condensate at finite 
temperature and found that for suitable initial velocities 
the impurity acquired energy, on average, after the first 
collision. However processes with multiple collisions were 
not considered. Here we show that the maximum cool- 
ing per impurity atom is obtained, for T not too small, 
at thermal equilibrium, that is, with multiple collisions. 

To simplify the problem we consider an homogeneous 
condensate, in fact it is meaningful in a trapped system 
in the sense of the Thomas-Fermi approximation. If the 
number of outcoupled atoms is sufficiently low, the initial 
thermal energy of the impurities can be neglected. The 
first cooling is obtained by the Raman scattering, that 
creates a cloud with a smaller temperature ||. We first 
introduce a Boltzmann equation for an impurity atom. 
A similar equation in momentum space is obtained, for 
a single particle and T = 0, in Ref. [?]. It is easy to 
show that the system goes toward the Boltzmann equi- 



librium distribution. We extend our discussion to a cloud 
of impurity atoms introducing suitable bosonic stimula- 
tion terms. We obtain an equation of the same type of 
the quantum Boltzmann equation |8|, ^|, [l0| . It is verified 
that the Bose-Einstein distribution is the steady state 
solution. 

Let us consider an impurity atom in an uniform con- 
densate at temperature T. If the impurity energy is 
Ei, the probability that the particle scatters in a state 
with energy E f is T{E f ,Ei) = T 1 {E f ,E l ) + Y 2 {E f ,E l ), 
where [|| 
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hLuQ(hujQ + 2/i) being, respectively, the initial impurity 

momentum, the condensate temperature in energy unit, 
the static structure factor and the energies of a free parti- 
cle and of a Bogoliubov quasiparticle of momentum q. M , 
fi and a are the atomic mass, the chemical potential and 
the scattering length for s-wave collisions between the im- 
purity atoms and the condensate ones, /i is equal to nog, 
where no is the condensate density and g = AitTi 2 ao / M , 
ao being the condensate s-wave scattering length. 

Fi and T 2 are associated with dissipative and cool- 
ing processes, respectively. In the first case condensate 
phonons are created, in the second one they are annihi- 
lated [p|. With a little of algebra we find that T(Ef, Ei) — 



2('y/fj,)G(2Ef/fi,2E i /iJ,), where 7 = ATm a c, c 
' n/M being the sound velocity, and 
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with ft = (3^/2 and Ae = ef — Ci. W(ef, ej) is 1 when > 
£/ > —1 when e/ > > 1/ef and zero elsewhere. 

The Boltzmann equation is 
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FIG. 1: Average energy acquired by the impurity in feT 
unit as a function of r, for some values of ei n . In the inset 
is reported the relative density Pi(r) as a function of r. (a) 
1/(3 = 10, i.e. k B T = 5/i. (b) 1/(3 = 2, i.e. k B T = /i. 



where P(2E/fi) is the density probability per unit energy 
and t = jt. It is easy to demonstrate that P(e) oc y/ee^^^ 
is a steady solution of Eq. (0) . 

The system dynamics depends exclusively on the prod- 
uct /3[x — ^/(kgT) and the initial probability distribu- 
tion. The condensate density and the scattering length 
are present in the rate factor 7 and, therefore, they in- 
fluence only the time scale of the dynamics. 

In Fig. 0a we report the average energy acquired by 
the impurity and lost by the condensate thermal atoms, 
in k B T unit, i.e. (3[(e) - e in ] = (3[{E) - E m ]. We have 
considered 1/(3 = 10 and some values (> 1) of the initial 
impurity energy. We can note that (E) = (3/2)fcsT for 
f — ► 00, that is the average energy of a free particle at 
the equilibrium. In contrast, if we neglect the collective 
modes, a thermal particle of the condensate has an av- 
erage energy equal to ~ 0.514(3/2)/csT, that is, about 
one half of the classical value. It is easy to demonstrate 
that the number of lost thermal atoms SN, when the 
condensate energy decreases by 5E, is ~ 0.78SE/(kBT). 
Therefore , if the impurity acquires an energy of ~ 1 . 3 k B T 
a thermal atom jumps into the condensate state, bal- 
ancing the loss due to the impurity creation. A nearly 
complete thermalization occurs for r ~ 0.5, that is, for 
tth — l/(2cronoc), where er = 47ra 2 . For a 23 Na con- 




densate and a typical density of 10 14 -7- 
350 + 3500^ for 1/(3 = 10) we have t th 



10 15 cm~ 3 (T - 
~ 10-=- 0.3ms. 



FIG. 2: Average energy acquired by the impurity in fi/2 unit 
as a function of r, for some values of tin and for 1/(5 — 0.6-rT. 



We can see that the thermalization time is independent 
on the initial velocity. In fact, the impurity loses quickly 
the memory of its initial energy. In the inset of Fig. 0a we 
report the relative density Pi(r) = P(ei n , r)/P(ej n , 0) as 
a function of r. The initial state population decreases by 
60 -j- 70 per cent for r ~ 0.1, i.e. for t — 0.06 2ms. We 
note that the depletion rate is smaller for a higher initial 
velocity, contrary to the zero temperature case. In fact, 
for T — we expect that the scattering rate grows upon 
increasing the initial velocity. At finite temperatures, 
to reduce the laser mode depletion, it is convenient to 
have a large initial velocity, because the scattering rate 
is smaller and, furthermore, the impurity leaves the con- 
densate earlier. For a sufficiently high velocity we expect 
that the scattered impurity fraction is not enhanced by 
the finite temperature. 

If the condensate temperature is reduced, the scatter- 
ing rate decreases and approaches the zero temperature 
value. For T = only the spontaneous scattering con- 
tributes and when ej„ < 1 the impurity thermalization 
is absent because of the condensate superfluidity. For 
ksT = \l {(3 = 2) the temperature can be still important. 
In fact, we can see in Fig. l]b that a nearly complete ther- 
malization occurs for r ~ 2, that, for hq — 10 15 cm~ 3 , 
corresponds to t ~ 1.2ms. The initial velocity is be- 
low the Landau value, so the first impurity scattering is 
due to the finite temperature. The function G(e/,ei) is 
equal to zero for £/ < 1/e;, therefore, if the initial en- 
ergy is too low, the particle jumps in a state with a large 
energy. However the probability rate of this process is 
small, because of the presence of the exponential in G. 
So we have a slowing down of the thermalization rate 
(see Fig. [l]b, solid line) when the initial velocity is very 
low. For ei n > 0.4 the cut-off of the constraint e/ > 1/e.; 
has a minor effect and the thermalization time becomes 
again independent on the initial velocity. Also the de- 
pletion rate of the relative density Pi(r) is reduced by 
the cut-off for ej n = 0.2. Increasing ej n the depletion 
rate grows, however for ej„ > 0.6 it slows down again, 
as it occurs for (3 = 10 (see Fig. |j]a). For ti n = 0.4 
and no = 10 14 -7- 10 15 cm -3 the initial impurity velocity 
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is ~ 0.3 -7- lcm/ 's, therefore, the average distance before 
the first collision (r ~ 1) is about 4 • (10~ 2 -j- 10~ 3 )mm. 
The radius of a condensate with the same central density 
is ~ 0.5 -j- 2 • 10 mm when the trap frequency is 20Hz. 
To contrast the gravitational acceleration of the impuri- 
ties in the radial direction of a horizontal cylindrical trap 
we can use a gaussian red-detuned laser beam. We can 
set this beam and the magnetic field in a such way to 
permit that the impurities leave the condensate after the 
thcrmalization time. 

In Fig. |2| we report the same curves of Fig. [ij but the 
energy is in p/2 unit and 1//3 = 0.6 1. In this case 
the energy reaches a maximum at an intermediate time. 
For 1/(3 — 1 and = 0.6 -f- 1.0 the maximum is reached 
at r = 10 -j- 20, that, for no = 10 15 cm~ 3 , corresponds 
to t = 6 -7- 10ms. We find that the evolution of the 
initial state population is nearly independent on ej n . At 
t = 10 the initial state is depleted by 40%. To cool the 
condensate it is not necessary to reach the maximum, 
but for smaller times we obtain a minor reduction of the 
condensate temperature per atom. Note that for 1//3 = 
0.6 the equilibrium thermal energy is 0.9, in p/2 unit. 
Therefore, for ej n = 0.9 4- 0.95 the initial energy is equal 
or higher than the equilibrium value. Even so, at an 
intermediate time the impurity acquires energy and the 
condensate cools. That occurs because when the velocity 
is below the critical value the first scattering can only 
enhance the impurity energy. 

Let us suppose now that there are many impurities, 
then many atoms can be in the same energy level and, 
therefore, it is necessary to consider the stimulation ef- 
fects. If we suppose that the initial impurity velocity 
is isotropic, we can make the following modification in 
Eq. (@) 
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where p(e) = P(e)/V, V being the system volume, and 
p{2E/n) = {2m) 3 / 2 h~ 3 (2ir)- 2 y/E is the level density per 
unit volume [jll| [l^]. It is easy to demonstrate that the 
Bose-Einstein distribution p{e) cx y / e(e^ e_ W — 1) _1 is a 
steady solution of Eq. (||) , p being the normalized chem- 
ical potential of the impurity cloud. 

We isolate in p(e) the contribution of the initial popu- 
lated state with energy ej n . The associated distribution, 
that is very narrow, has a width equal to Se = 2/(Vpp). 
We indicate with its height. 

Eq. (H) becomes 
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FIG. 3: Relative density Pi(r) as a function of r for 1//3 = 2 
and some values of ni/n t h, with (a) ei n = 0.2 4 0.4, (b) = 
1.5 and e™ = 1 in the inset. 



The quantities SeG(ei n , e)p(e) and J G(ei n ,e)p(e)de in the 
first and second equation, respectively, can be neglected. 
Therefore Eqs. become 

^ = B{p) + -{\p(e)/p(e) + 1]G(6, e m ) 
dr p 

-G(e in ,e)p{e)/ p(e in )}ni (6) 
= -7t(p)»<) ( 7 ) 
where is the spatial density of the impurity cloud and 

liip) = I \p(e)/p(e) + l]G(e, e in )de - 
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l/p{e in ) / G(e in ,e)p(e)de = 



p(e) de,(8) 



y/1 + Ae 2 /4 



Ae being Ci„ — e. The first term of 7, is the depletion rate 
of the initial state when the stimulated effects of the im- 
purities are negligible. The second term, that does not 
depend directly on the temperature, is proportional to 
the density p(e) and it is important when there are many 
impurities. This term is negative when ei„ < 1, therefore 
in this case it reduces the depletion rate 7^. In fact, the 
scattered atoms can have a high probability to scatter 
back into the original macroscopically populated state. 
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FIG. 4: Same as with Fig. |lj but for an impurity cloud. Some 
values of m are considered, for 1/(3 = 2 and e in = 0.6. 



On the contrary, when ei n > 1 the second integrand is 
negative for e > ej„ and positive for e < ej n ; if ej„ is suffi- 
ciently high the impurity stimulation effect enhances the 
depletion rate, as it has been experimentally observed 0. 
In our calculations we have considered the condensate as 
a thermal reservoir, however if we have many impurities 
the collisions can cool or warm the condensate, conse- 
quently the depletion rate can decrease or increase. Our 
approximation is reasonable when the density of outcou- 
pled atoms is small with respect to the density of thermal 
atoms n t h or, otherwise, for sufficiently small times, when 
the collided fraction is small. We have approximatively, 
if we neglect the condensate collective modes, Jll|] 



nth = (4 7 r^ 3 )- 1 m 3 / 2 ( M //3) 3 / 2 



(9) 



In Fig. ||a we report the relative density P%{t) for = 
2, e in = 0.2-^0.4 and some values of nj/nth = n>i(0) /nth- 
It is easy to demonstrate that these curves are indepen- 
dent on /i and, therefore, on a and uq. In fact, with 
p(e) cx JJi and rij cx fJ 3 ^ 2 , the Boltzmann equations are 
invariant with respect to modifications of the chemical 
potential. For no = 10 14 cm~ 3 and uq 



have nth — 1.1 • 10 



and n t h — 3.5 • 10 13 cm 3 , re- 



spectively. We can see that for ni/n t h — 10 the depletion 
rate is considerably reduced by the stimulation effects. In 



Fig. ||b we plot the same curves but for e,„ = 1 (inset) 
and 1.5, respectively. These curves are reliable only at 
the initial times because of the cooling of the condensate. 
Furthermore, if collisions create many thermal impurity 
atoms, a lot of impurities could condensate into the state 
with zero energy, but our equations are not suitable to 
describe this phenomenon, because this state has to be 
treated separately from the energy integral of Eq. (||). 
If the number of scattered atoms is lower than n t h we 
can suppose that this condensation does not occur. Note 
that the depletion rate is lowered also for ej„ = 1.5, that 
is above the Landau value. Instead, for ej„ = 2 we find 
that the depletion rate is enhanced. Increasing the tem- 
perature the threshold of €i n increases. 

Finally, we report in Fig. |] the same curves of Fig. [I], 
but accounting for the stimulation effects of the impurity 
cloud. We have considered some values of rir, for 1//? = 2 
and ei n = 0.6. The maximum acquired energy is reduced 
increasing nj. 

In conclusion, we have evaluated the thcrmalization 
time of impurity atoms inside a condensate at finite tem- 
perature and we have found that their thcrmalization, 
that occurs also below the Landau critical velocity, can 
cool the condensate. In fact, a small Raman or rf outcou- 
pling creates a condensate in another level with a lower 
temperature ||, that cools sympathetically the parent 
condensate. This effect can be studied experimentally 
creating untrapped impurities whose gravitational accel- 
eration is compensated by a detuned laser field. The cool- 
ing effect could be used to remove the residual thermal 
atoms. Furthermore, we have accounted for the stimula- 
tion effects of the impurity cloud and we find that below 
the Landau critical velocity the macroscopic population 
of the initial state can reduce drastically its depletion 
rate. For sufficiently high velocities the opposite effect 
occurs, as observed in Ref. |2[ The threshold velocity 
depends on the condensate temperature. Our analysis 
is also important for the development of intense atom 
lasers. 
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